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Abstract 

We study the combined effect of temperature and disorder in the spacer on the interlayer 
exchange coupling. The temperature dependence is treated on ab initio level. We employ 
the spin-polarized surface Green function technique within the tight-binding linear muffin-tin 
orbital method and the Lloyd formulation of the IEC. The integrals involving the Fermi-Dirac 
distribution are calculated using an efficient method based on representation of integrands 
by a sum of complex exponentials. Application is made to Co/Cuioo-xM x /Co(001) trilayers 
(M=Zn, Au, and Ni) with varying thicknesses of the spacer. 

1 INTRODUCTION 

The oscillatory interlayer exchange coupling (IEC) between magnetic layers separated by a non- 
magnetic spacer has recently attracted considerable attention in the literature. The physical 
origin of such oscillations is attributed to quantum interferences due to spin-dependent confine- 
ment of electrons in the spacer. The periods of the oscillations with respect to the spacer thickness 
are determined by the spacer Fermi surface, and this conclusion has been confirmed by numer- 
ous experiments. In particular, the change of the Fermi surface by alloying thus leads to the 
change of the oscillatory periods (van Schilfgaarde et al, 1995; Kudrnovsky et al, 1996). On the 
other hand, there are very few studies of the temperature-dependence of the IEC (Bruno, 1995; 
d' Albuquerque e Castro et al, 1996), and its systematic study on ab initio level is missing. 

The main mechanism of the temperature dependence of the IEC is connected with thermal 
excitations of electron-hole pairs across the Fermi level as described by the Fermi-Dirac function. 
It turns out that other mechanisms (e.g. electron-phonon or electron-magnon interactions) are 
less important. The effect of the temperature on the IEC can be evaluated either analytically or 
numerically. The analytical approach assumes the limit of large spacer thicknesses, for which all 
the oscillatory contributions to the energy integral cancel out with exception of those at the Fermi 
energy. The energy integral is then evaluated by a standard saddle-point method (Bruno, 1995). 
The general functional form of the temperature-dependence of the interlayer exchange coupling 
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S X (T) in the limit of a single period is then: 

cNT 

£x(T) = £ x (0)t(N,T) , t(N,T) = — j^Ff) • ID 

Here, T denotes the temperature, A" is the spacer thickness in monolayers, and c is the con- 
stant which depends on the spacer Fermi surface. The term S x (0) exhibits a standard N~ 2 - 
dependence (Bruno, 1995), while the scaling temperature factor t(N,T) depends on N via NT. 
In the preasymptotic regime (small spacer thicknesses) the functional form of t(N, T) differs from 
that of Eq. (TO), particularly in the case of the complete but relatively weak confinement due to the 
rapid variation of the phase of the integrand which enters the evaluation of the IEC (d 'Albuquerque 
e Castro et ai, 1996). 

The second, numerical approach is in principle exact, not limited to large spacer thicknesses, 
however, it may be numerically very demanding, in particular for low temperatures. It is applicable 
also to disordered systems with randomness in the spacer, magnetic layers, or at interfaces (Bruno 
et ai, 1996). 



2 FORMALISM 

The multilayer system consists of the left and right magnetic subspaces separated by a non- 
magnetic spacer (the trilayer). The spacer may be a random substitutional alloy. We employ the 
Lloyd formulation of the IEC combined with a spin-polarized surface Green function technique as 
based on the tight-binding linear muffin-tin orbital (TB-LMTO) method. The exchange coupling 
energy £ X (T) can be written as 

£ X (T) = Im/(T), I{T)= [ f{T,z)*(z)dz, (2) 

where f(T, z) is the Fermi-Dirac distribution function and 

= J^E^ lnM ( k ll^) (3) 

is a difference of (in the case of disorder, of configurationally averaged) grandcanonical potentials 
for the antiferromagnetic and ferromagnetic alignments of magnetic slabs (Drchal et ai, 1996). 

The energy integration is performed over a contour C along the real axis and closed by a 
large semicircle in the upper half of the complex energy plane, tr^ denotes the trace over angular 
momentum indices L = (£m), the sum runs over ky-vectors in the surface Brillouin zone, and Nu is 
the number of lattice sites in one layer. The quantity M (kii , z) is expressed in terms of the screened 
structure constants which couple neighboring (principal) layers and of the so-called surface Green 
functions. All details can be found in (Drchal et ai, 1996). We only note that the use of a Green 
function formulation of the IEC is essential for describing the randomness in the spacer within 
the coherent potential approximation (CPA) which is known to reproduce compositional trends 
in random alloys reliably. 

The integral in (0) can be recast into a more suitable form using the analytic properties of 
^f(z), namely, (i) ^f(z) is holomorphic in the upper complex halfplane, and (ii) z^f(z) — > for 
z — > oo, Imz > 0. Let us define a new function $(?/) = — i ^/(E F + iy) of a real variable y, y > 0. 
Then at T = K, 

r+oo 

1(0) = $(y)dy, (4) 

JO 
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while at T > K, 



I(T) = 2nk B T y £ l $(y k ), (5) 



k=l 



where k B is the Boltzmann constant and y k are Matsubara energies y k = nk B T(2k — 1). In the 
limit T — > 0, I(T) — > 1(0) continuously. 

We have verified that the function $(y) can be represented with a high accuracy as a sum of 
a few complex exponentials in the form 

M 

Hv) = E Ajexpfey), (6) 
i=l 

where Aj are complex amplitudes and pj are complex wave numbers. An efficient method of 
finding the parameters Aj and pj is described elsewhere (Drchal et al, 1998). The evaluation of 
I(T) is then straightforward: 

M A . 

J(T) = -2ttA; b T £ J , (7) 

exp (nk B Tpj) - exp [-irk B Tpj) 



which for T = K gives 



i=i ^ 



3 RESULTS AND DISCUSSION 

Numerical studies were performed for an ideal fcc(OOl) layer stack of the spacer (Cu) and magnetic 
(Co) layers with the experimental lattice spacing of fee Cu. The spacer layers can contain the 
impurities (Zn, Ni, and Au) which form the substitutional alloy with spacer atoms. Possible lattice 
and layer relaxations are neglected. Alloying with Ni, Zn, or Au alters the electron concentration 
and, consequently, modifies the Fermi surface and thus, in turn, also the temperature dependence 
of the IEC. The most obvious effect of the alloying, for T = 0, is the change of the periods of 
oscillations connected with the change of the corresponding spanning vectors of the alloy Fermi 
surface (Kudrnovsky et al, 1996). The more subtle effect of the alloying is connected with damping 
of electron states and relaxation of symmetry rules due to alloying. 

To determine the parameters of complex exponentials @, we have evaluated the function 
$(?/) at 40 Matsubara energies corresponding to T = 25 K. We have verified that the results 
depend weakly on the actual value of the parameter T. Special care was devoted to the Brillouin 
zone integration. The efficiency of the present approach allows us to perform calculations with 
a large number of ky-points in the irreducible part of the surface Brillouin zone (ISBZ). Note 
also that such calculations have to be done only once and then the evaluation of the IEC for any 
temperature is an easy task. In particular, we employ typically 40000 ky -points in the ISBZ for the 
first Matsubara energies close to the Fermi energy and the number of ky -points then progressively 
decreases for points distant from the real axis. The present calculations agree with the results of 
conventional calculations (Drchal et al., 1996) but they are much more efficient numerically, in 
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particular when calculations for many different temperatures are required. 




Figure f: Absolute values of the discrete Fourier transformations of N 2 £ X (N,T) with respect to 
the spacer thickness iV as a function of the temperature T and the wave vector q for a trilayer 
consisting of semiinfinite Co-slabs sandwiching the spacer with indicated compositions. 

The calculations were done for spacer thicknesses N — 1 — 50 monolayers and for temperatures 
T = — 500 K (in steps 10 K) and by assuming semiinfinite Co-slabs. In this case only one period, 
namely the so-called short-period exists, which simplifies the study. There are several possibilities 
of how to present results (see Drchal et al, 1998) for more details). As an illustration, in Fig. 1 
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we plot the discrete Fourier transformations (Drchal et al., 1996) of N 2 S X (N,T) with respect to 
N, S x (q,T), as a function of variables q and T. The discrete Fourier transformation on a subset 
N G 10 — 50 which avoids the preasymptotic region is employed here. The positions of peaks 
of q = q m then determine oscillation periods p = 2n/q m , while \£ x \ give oscillation amplitudes 
{Drchal et al., 1996). In particular, one can see how the modification of the Fermi surface due to 
alloying changes the temperature dependence of the IEC, i.e., the coefficient c in Eq.(|TJ). 

The following conclusions can be drawn from numerical results: (i) The non-random case (Cu) 
exhibits the period p ~ 2.53 MLs (monolayers) or, equivalently, q m pa 2.48 in accordance with 
previous calculations (Drchal et al, 1996). In accordance with (Kudrnovsky et al, 1996), the 
periods of oscillations for Cu75Zn 25 alloy are shifted towards higher periods (p ~ 3.05MLs), to- 
wards smaller periods for CU85NI15 alloy (p ~ 2.27 'MLs), while they remain almost unchanged 
for equiconcentration CuAu alloy spacer (p pa 2.36 MLs); (ii) The periods of oscillations are 
temperature independent because the electronic structure or, alternatively, spanning vectors are 
temperature independent; (iii) The amplitudes exhibit a strong temperature dependence in agree- 
ment with predictions of model theories (Bruno, 1995). In particular, our results agree reasonably 
well with those of Fig. 3 in (d 'Albuquerque e Castro et al., 1996) for the case of ideal Cu spacer, 
(iv) For alloy spacers at T = we mention, in particular, the dependence N~ 2 of the oscillation 
amplitudes on the spacer thickness N for CuNi and CuAu alloy spacers, while additional expo- 
nential damping due to disorder was found for the CuZn alloy. This indicates a finite lifetime 
of states at the Fermi energy for ky-vectors corresponding to short-period oscillations for this 
case and only a weak damping for CuAu and CuNi alloys; (v) Finally, the effect of temperature 
(the factor t(N,T) in Eq. (|l])) is similar for a pure Cu-spacer, CuNi and CuAu alloys, but it is 
much smaller for the case of CuZn alloys. The effect of temperature, similarly to alloying, is to 
broaden spanning vectors of the Fermi surface (Bruno, 1995). If the damping due to alloying 
is non-negligible, the combined effect of disorder and temperature leads to a relatively smaller 
(compared to the case T =0 K) suppression of the oscillation amplitude with the temperature. 
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